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ABSTRACT 

We present joint constraints on the distribution of Mg II absorption around high redshift galaxies that are 
obtained by combining two orthogonal probes, the integrated Mg II absorption seen in stacked background 
galaxy spectra and the distribution of parent galaxies of individual strong Mg II systems as seen in the spec- 
tra of background quasars. We present a suite of models that can be used to predict, for different two and 
three-dimensional distributions, how the projected Mg II absorption will depend on a galaxy's apparent incli- 
nation, the impact parameter b and the azimuthal angle between the projected vector to the line of sight and the 
projected minor axis. The models are all constrained to have the same average projected radial profile, charac- 
terized by a power-law with an exponential cut-off at Rcut = 37.5 kpc. In general, we find that variations in 
the absorption strength with azimuthal angles provide much stronger constraints on the intrinsic geometry of 
the Mg II absorption than the dependence on the inclination of the galaxies. In addition to the clear azimuthal 
dependence in the integrated Mg II absorption that we reported earlier in Bordoloi et al (2011), we show that 
strong equivalent width Mg II absorbers (Wr(2196) > 0.3 A) are also asymmetrically distributed in azimuth 
around their host galaxies: 72% of the absorbers in Kacprzak et al (2011), and 100% of the close-in absorbers 
within 38 kpc of the center of their host galaxies, are located within 50° of the host galaxy's projected semi 
minor axis. It is shown that either composite models consisting of a simple bipolar component plus a spherical 
or disk component, or a single highly softened bipolar distribution, can well represent the azimuthal dependen- 
cies observed in both the stacked spectrum and quasar line datasets within 40 kpc. Simultaneously fitting both 
data sets, we find that in the composite model the bipolar cone has an opening angle of ~ 100° (i.e. confined 
to within 50° of the disk axis) and contains about 2/3 of the total Mg II absorption in the system. The single 
softened cone model has an exponential fall off with azimuthal angle with an exponential scale length in open- 
ing angle of about 45°. We conclude that the distribution of Mg II gas at low impact parameters is not the same 
as that found at high impact parameters. Mg II absorption within 40 kpc primarily arises from cool Mg II gas 
entrained in star formation driven winds. At larger impact parameters beyond 40 kpc, there is evidence for a 
much more symmetric distribution, significantly different from that closer in to the galaxies. At larger radii, a 
significant component appears which is aligned more with the disk and is possibly inflowing, perhaps as part 
of a galactic fountain or the inflow of material from further out in the system. 

Subject headings: galaxies: evolution: general — galaxies: high-redshift — intergalactic medium — ISM: jets and 
outflows — quasars: absorption lines 



1. INTRODUCTION 

The absorption line systems observed in the spectra of 
background sources such as quasars and galaxies are tracers 
of the otherwise invisible gaseous structures around galaxies. 
Combining such absorption-line observations with informa- 
tion on their host galaxies, provide a unique and powerful 
means to understand the baryonic processes in and around 
galaxies. 

The Mg II AA 2796, 2803 doublet is one of the most stud- 
ied tracers of the gaseous environment around galaxies. These 
absorption lines are believed to originate in p hoto ionized gas 



at temperatures of around T^ 10"^ K ( [Bergeron & Stasinska 
|1986t Charlton et al.|[20Q3] ) and to trace a wide range of 



neutral hydrogen column densities of Nhi ^ 10^^- 10^^ cm ^ 
( [Churchill et aLl[2QQQl [Rigby et al.|[2Q^ [Rao et al.|[2006l ) 
within a few hundred kpc of their host galaxies. 
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There have been many detailed studies on the distribution 
of column densities, the redshift evolution of number densities 
and the kinematic signatures of Mg II absorbers using quasar 
absorption line systems (see e.g. Lanzetta et al. 1987; Sar 

Sar- 



gent et al.|[T 988*, 'Pe titjean & B ergeron 1990> ,Steidel & Sar- 
gent[[1^^2l ^C harlton & (ljhil? 8hill 1998[ [Nestor et ai.[|^00^[ 
^dT • • ■ 



Prochter et al.n2006| . The association of strong Mg II ab- 
sorption with normal, bright, field galaxies is by now well 
established (e.g. Churchi ll et al.|20 05 and references herein). 
The Mg II gas around galaxies is traced out to ^ 100 kpc with 
absorber covering fractions of 50-80% ( Ka cprzak et al. 12008* 



[Chen et al. [20 lOj [Nielsen et al.[2012 ). The anti-correlation be 
tween absorber equivalent width and impac t parameter is also 
now well established by several studies ( Steideri995VBouche 
et al. 2006, , Chen et al. 2010, Kacprzak et al. 201 lb, Bordoloi 



etar.|201l| r 



While the association of Mg II absorption with galactic 
haloes has been secure for some time, the origin and fate 
of this gas has been less clear. Possibilities have included 
outflowing material entrained in star-formation driven winds 
([Weiner et al. 2009; Rubin et al. 2010 ) or inflowing material 



feeding the disks of galaxies (Rubin et al.,201 2. Martin et al. 
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The process of large scale outflows driving cold gas out of 
the galaxy and polluting the circum-galactic medium (CGM) 
is certainly a complex one. Moreover, the kinematics and the 
distribution of cold gas entrained i n such s tar formation driven 
outflows still remain uncertain ( Veilleux et al. 2005 ). The 
mechanical energy generated by supernovae explosions and 
cosmic rays etc. create hot over-pressurized bubbles which 
sweep up the surrounding ambient material into dense shell 
like structures. These shells accelerate as they expand and 
when the reach the low density CGM, Rayleigh-Taylor insta- 
bilities set in, which leads to fragmentation of these shells in 
the CGM (Heckman 2002; Fujita et a l. 2009 ). Although the 
physical processes involved in such a scenario are likely to 
be complex, it is likely that cold pockets of gas will be en- 
trained in such outflowing material. Moreover, new clumps 
of cold gas m ight form due to thermal instabilities in the hot 
gas ( Fujita et al. 2009 ). Of course, entrained material ca n also 
later fall back onto the galaxy due to gravity ( Oppenheimerl 
et al.|2010 ). Moreover, derived constraints on mass, energet- 
ics or momentum of the outflow are hindered with orders of 
magnitude uncertainties ([Prochaska et al.|20TT]and references 



see Churchill et al. ( 2012 ), who find the same result using 51 
galaxies with sensitivities to Wr(2796) > 0.003 A. 



therein). Even idealized outflow models (Fu jita et al.J|2009| ) 
are not well constrained. The opening angle of the outflowing 
gas, the acceleration of the gas as a function of galactocentric 
radius, optical depth, covering fraction, density, temperature 
or the spatial extent of such outflowing gas are all quite uncer- 
tain. For more complex models, where the Mg II absorption 
occurs due to combination of outflowing gas with extended 
galactic disk or inflowing gas, the fractional contribution from 
each component also remain unconstrained. 

Recently, several studies have sought to infer the distribu- 
tion of Mg II absorption around galaxies so as to provide diag- 
nostic information on the origin of the enriched material and 
to try to constrain possible models. Using very large numbers 
of about 5000 ba ckground galaxy spectra from zCOSMOS 
( [Lilly et aL||2007| ) to probe the profiles of Mg II a bsorption 



around 4000 foreground galaxies at0.5<z<0.9, [Bordoloi 
[et al.[ ( [201 1| ) showed that the strength of the Mg II absorp- 
tion was much higher for blue star-forming galaxies than for 
red, more passive, galaxies, at the same stellar mass. Further- 
more, they showed that for inclined foreground disk galax- 
ies, the Mg II absorption is significantly enhanced along the 
galaxy projected minor axis, at least out to 40-50 kpc from the 
galaxy. This result indicated that the Mg II absorption must 
be primarily originating in a bipolar structure aligned with the 
disk axis, suggesting that it is part of a bipolar outflow driven 
by star formation. 

On the other hand, [Kacprzak et al.| ( [201 Ibj ) found an anti- 
correlation between Mg II absorber strength and disk inclina- 
tion /, as normalized by the impact parameter h, in the sense 
that the equivalent width of Mg II was anti-correlated with 
b/i. This was interpreted as being due to co-planner geom- 
etry of the absorbers. In parallel Bouche et al. ( 2011), using 
the sample of ten z ~ 0. 1 absorbers of Kacprzak et al. ( 201 la) 
found a bimodal distribution of azimuthal angles, 60% of 
their absorbers being within 30° of the galaxy minor axis. 
They also interpret their re s ult as to support the bipolar wind 



scenario. 



Kacprzak et aL] ( [2012 ), using 88 1^,(2796) > 0.1 
A absorbers and 35 Wr(2796) < 0.1 A non-absorbers also 
reported a bimodal distribution of azimuthal angles: They 
found that the equivalent width distribution for Mg II absorp- 
tion along the disk major axis are likely skewed towards the 
weak absorbers than those along the disk minor axis. (Also 



The observational material present ed in Bordoloi et al. 
([201 1[) and in the studies of Kacprz ak et al.[ ( [201 lb| ) ana 
Bouche et al. (2011) are in many ways orthogonal. The in- 
tegrated approach using stacked galaxy spectra measures the 
total Mg II absorption, integrated over all equivalent widths, 
in some location relative to the host galaxies. This location 
can be defined over a range of impact parameters h and/or az- 
imuthal angles (which we throughout the paper define to be 
relative to the projected minor axis). By stacking the spectra 
to detect the absorption signal in the first place, this total Mg 
II absorption is by construction measured as an average over 
some potentially large sample of foreground galaxies that can 
be constructed as desired in terms of their luminosities, col- 
ors, masses, star-formation rates, apparent inclinations and so 
on. 

In contrast, the quasar line approach involves identifying 
absorption systems above some individual equivalent width 
threshold, and then identifying the host galaxies and deter- 
mining both the location of the quasar line of sight relative to 
them, and any properties of the host galaxies of interest (col- 
ors, star- formation rates etc.). The distribution of these mea- 
surements relative to what would have been expected from a 
sample of galaxies without the Mg II selection can then give 
information on the Mg II itself. A specific example will make 
the distinction clear. In the stacked approach, one measures 
directly the total average Mg II absorption strength as a f{(j)) 
around galaxies. In the quasar-line approach, one measures 
the observed distribution of A^(^) for the hosts of detected Mg 
II systems and compares this with what would be expected 
from the general population, i.e. in this case, a flat distribu- 
tion in N{(l)). 

The main aim of the current paper is to critically examine 
the information that comes from both the integrated spectra 
approach and the quasar-line approach to see if they present a 
consistent picture for the distribution of Mg II around galax- 
ies. Then, since we find that they are indeed consistent, we 
wish to undertake a joint analysis, using both data sets, to 
define, as well as possible, the global characteristics of this 
distribution. For both purposes, we need a clear understand- 
ing of the role of the inclination of a galaxy in determining the 
Mg II absorption as seen projected onto the plane of the sky. 
We therefore construct an extensive suite of models that can 
be viewed from any angle and with any impact parameter. 

The layout of this paper is as follows. In Section 2 we first 
discuss the construction of simple geometric models for Mg 
II around galaxies, including the computation of an integrated 
equivalent width and the definition of the spatial distributions 
and geometries of the different models that are considered for 
the simulations. In Section 3 we "observe" these models over 
the full range of impact parameter b, inclination /, and az- 
imuthal angle (/), both to identify the most powerful diagnos- 
tics for the geometry and to provide predictions for the run 
of absorption strength with these observable quantities that 
can be used to compare with the observations. In Section 
3.2 we compare these model predictions to the stacked in- 
tegrated Mg II measurements, including considering straight- 
forward superposition of different components. In Section 4 
the model predictions are compared with the quasar absorp- 
tion line data, yielding similar results. In section 5, we per- 
form a joint analysis of both the datasets to derive our globally 
preferred model. Our findings are summarized in Section 6. 
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Fig. 1 . — Top panel: Cross section of the Nui density distribution (in the plane of the sky) for the sharp bipolar cone (top left) and for the softened bipolar cone 
model (top right) both with an opening angle of 45° viewed edge on. Bottom panel: The projected Nhi column density distribution for a 2D disk viewed at an 
inclination angle. For all three cases, the models have the same projected radial profile, with Rent - 37.5 kpc. 



Throughout this paper, we adopt a ACDM cosmology with 
Vim = 0.27, l^A = 0.73 and = 70 km s'^ Mpc'K 

2. THE MODEL 

We will assume throughout the paper that the light from the 
galaxy is in a disk and that the ellipticity of this light indicates 
the inclination of the disk. We describe this orientation by 
the direction of the pole of this disk, which lies in the direc- 
tion P. We then create a set of galaxy models that are each 
based on a particular 3 -dimensional distribution for the ab- 
sorbing gas that is given in terms of a density niR) which is 
the number density of the neutral hydrogen gas (number of 
atom per cnT^) in the galaxy at position R. The angle sub- 
tended by R with respect to the polar axis P is defined to be V^, 
i.e. cos = R.P. We consider a spherically symmetric distri- 
bution, a bipolar distribution and a two-dimensional disk. In 
each case the density distribution is assumed to be rotationally 
symmetric about the pole. 

The galaxy is then rotated randomly through 10^ orienta- 
tions, i.e. so that P is uniformly distributed on the unit sphere. 



The expected Mg II absorption is then calculated by integrat- 
ing through a given line of sight that passes at a projected sep- 
aration, i.e. impact parameter, b from the centre of the galaxy. 
The observed column density Nniib) at impact parameter b is 
given as 



NHi(b) 



= J n 

LOS 



(R) dl 



(1) 



Each of the 10^ random orientations is associated with a 
particular inclination for the observer, where the inclination 
angle / is defined so that / = 0° corresponds to a face-on 
galaxy and / = 90° corresponds to an edge-on galaxy whose 
pole lies in the plane of the sky. The inclination angle / has 
a probability distribution P(i)di ~ sin(i)di, i.e. more galaxies 
are edge-on than face-on. For simplicity, the lines of sight 
are chosen to all lie along a particular line in the sky extend- 
ing from the center of the galaxy, so that each orientation also 
corresponds to a projected azimuthal angle (j) between the pro- 
jected pole of the galaxy, i.e. the observed minor axis, and 
the arbitrary vector in the plane of the sky along which the 
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Fig. 2. — The azimuthally averaged radial profile for all models studied in 
this paper, which are by construction the same. For all models the overall 
normalization (a) is so c hosen that (W^^y) .7 A within 20 < ^ < 50 kpc. 
The data points are from |Bordoloi et al!] ( |2011) . 

lines of sight have been calculated. We use ^ = 0° for a 
sight line that is located along the minor axis, i.e. above the 
pole of the galaxy, and (j) = 90° corresponds to a sightline 
aligned along the major axis. This azimuthal angle (cj)) will 
be uniformly distributed between 0° < < 90° . Through- 
out this paper, kinematics of the Mg II absorbers are not taken 
in to account. We quote the total equivalent width of the ab- 
sorption absorption line system and assume no knowledge of 
their kinematic velocity spread. All the models wi ll be nor- 
maliz e d to th e radial profile of blue galaxies from Bordoloi 
1 ), hence the models studied here will represent 
blue host galaxies within the stellar mass range of ~ lO^^M© 
to IO^^Mq. 

For comparison with other works, we ran the models using 
neutral hydrogen number densities, n(HI) and column densi- 
ties Nhi- For a given galaxy orientation, we then convert the 
integrated column density Nniib) to an observed equivalent 
width in Mg II, Wobsib). We ta ke two approaches. First, we 
adopt the linear expression from|Menard & Chelouche|(|2009|) 



NHi{Wobs) = {2A5±03^)x\0 



19 



0.5 -W, 



obs 



2.08±0.24 



(2) 

which has no effects of line saturation. It should be noted that 
the Wobs values quoted here are summed over both the compo- 
nents of the Mg II absorption doublet (Wobs = ^2796 + ^2803)- 
Since we normalize our models to an observed radial profile 
Wobs(i^) the use of the HI number density and the adoption of 
the above conversion factor are simply book-keeping devices 
which have no impact on our results. 

Second, we modify this relation to account for line satu- 
ration. The effects of line saturation in the context of our 
model may be quite complex, depending on the cloud struc- 
ture of the Mg II absorbing gas. For instance, it would be 
possible to have saturated discrete clouds (producing satu- 
rated line ratios) whose integrated absorption is nevertheless 
added linearly along the line of sight through the galaxy be- 
cause the clouds have sufficiently different velocities. The 
much broader "beam" of studies using background galaxies 



rather than quasars is also relevant because small clouds could 
cover the pc^ area of the quasar beam but not the 10^ pc^ 
beam probed by a background galaxy. For these reasons, we 
adopt a simple ad hoc prescription for line saturation. Mg 
II absorption lines become saturated at column densities of 
^Mgii ^ 10~^^ cm~^, which is at relatively low hydrogen col- 
umn densities of Nhi > 10^^ cm~^ at solar abundance. There- 
fore we assume, for all lines of sight, that if Wobs(i^^ h (l))>2A 
as given by equation [2j then that line of sight is probing a sys- 
tem where the effects of saturation are likely to be present. A 
new modified Wobs,s is obtained by drawing from a Gaussian 
distribution with mean Wobs = 2.0 A and standard deviation 
cr = 0.3 A. The values are chosen to more or less approximate 
the scatter in the data used to estimate equation p] The ex- 
treme values of the new saturated Wobs^s, which fallbelow 0.6 
A are redrawn from a uniform random distribution between 
0.6 and 2.0 A. This is chosen so that the saturated systems do 
not correspond to unnaturally low equivalent widths. 

Since the quoted Wobs,s values are summed over both com- 
ponents of the Mg II 2796, 2803 doublet, we have to assume 
a doublet ratio to convert Wobs.s to the equivalent of Wr(2196) 
values. The oscillator strengths of the Mg II absorption dou- 
blet lines have a ratio of 2:1 and in the optically thin regime, 
their equivalent widths will also have a doublet ratio of 2:1. 
But as the optical depth of the lines increase and the lines 
become saturated, the equivalent width ratio also approaches 
1:1. In observed quasar absorption line systems the doublet 
ratios indeed vary from 2 to 1 for W^(2796) < 0.1 A. Further- 
more, doublet ratio is also a function of impact parameter. In 
this paper, we only consider strong (^^(2796) > 0.3 A) sys- 
tems, which will all be saturated. Hence for simplicity we 
assume a doublet ratio of 1:1, i.e. the Wr(2796) component 
will have half the value of Wobs- 

Given the radial profile {Wobs)(f^) for each of the 10^ ori- 
entations, each with their unique (/, (j)) we can then average 
these to yield observable quantities including (1) radial pro- 
files {Wobs)(i^), i-e. averaged over / & (/), (2) Wo/?^(/,0) at a 
given b, (3) the inclination profile : {Wobs)(i) for a given b, av- 
eraged over (j) and (4) azimuthal profiles: (Wobs) at given 
b and averaged over /. This exercise is then repeated for the 
different geometrical models constructed as described below. 

2.1. 3-D bi-cone and sphere models 

The first model is set up by assuming that absorbing gas 
is distributed with a conical geometry, aligned along the disk 
rotation axis. Within the galaxy, the angle with respect to the 
polar axis is defined to be ijj, i.e. cos?/^ = R.P. Within this 
cone, the number density is given by the galactocentric radius, 
R, and to be zero outside the cone. The radial dependence is 
set to be a power-law with an exponential cut-off, s ince this 
will give a good match to the average radial profile in jBordoloi 
et al..(2011 ). 



n(R) = { 



\ (£;)^exp(-£), ifV^<x/2 



(3) 



0, 



otherwise 



Here x is the opening angle of the bipolar cone. The generic 
models defined by equation|3]are hereafter referred to as sharp 
bipolar cone models. This model represents the simple model 
of bi-conical outflow coming out of a galaxy. We create a suite 
of models with opening angles x ^ [0°,180°]. The model 



Modeling the distribution of Mg II absorbers around galaxies 



5 




30 40 50 60 
Inclination Angle 




30 40 50 60 
Inclination Angle 




I 50- 

I 

I 40- 

N 
< 



20 30 40 50 60 
Inclination Angle 



Disk 




10 20 30 40 50 60 70 80 
Inclination Angle 




(0 

Q _ 

O) 

c 
o 

< - 



o 

Q. 

O) . 

c 
o 

< 



30 40 50 60 
Inclination Angle 



Face on along 
disk plane 



Highly inclined along 
disk plane 



All values at 
b = 45 kpc 



Face on along 
minor axis 



Highly inclied along 
minor axis 



' ' ' I ' " ' I 
Face on 



I ' " ' I ' ' ' ' I 
Edge on 



Fig. 3. — The \ogiQ{NHi(b)) column density distributions for different models as a function of inclination and azimuthal angle, for b = 45 kpc. The models 
shown here are respectively from top left (i) cone with opening angles (x) 45°, 90°, 130°, (ii) sphere and (iii) Disk. Different 3-d absorption geometries show 
clear trends with inclination and azimuthal angles. The last panel shows the orientation of the absorbers. 



with X = 180° is therefore a full sphere. 

We take (3 = -2.0 since this will w ell match the average ra- 
dial profile in [Bordoloi et al.| (201 shown also in figure |3] 



It is obviously appropriate for a mass-conserving outflow of 
constant velocity. The default cut-off radius in the exponen- 
tial (Rcut) is set at 37.5 kpc, also to match the average radial 
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profile. The value of a is set for each of the conical models 
so that the average {Wobs) ^ 0.7 A for 20 < /? < 50 kpc. This 
again is chosen to match the observed average radial profile in 
figure [3] Because the conical models all have the same radial 
dependence, there is a simple relation between the adopted a 
and the opening angle x- 

Top left panel of Figure [T] shows a cross section of the Nm 
column density distribution of the sharp bipolar cone model 
with X = 45° viewed edge on with / = 0. For the remainder of 
the paper we refer to sharp bipolar cone models with a sub- 
script indicating the opening angle x, i.e. M30 has an opening 
angle of 30° and Migo is the sphere. 

The second set of bipolar cone models are defined to have 
the same radial fall-off as previously described, but now with 
a softened exponential falloff in ?/^, i.e. 



R 

Rrul 



exp 



R 



Rr 



exp 



2^1) 
X 



(4) 



The generic set of models defined by equation |4] are here- 
after referred to as softened bipolar cone models. Again as 
before, we set Rcut = 37.5kpc and (3 = -2.0. The free pa- 
rameter a is set as before. This model has the feature of a 
non-zero density down to the plane of the galaxy and might 
represents a situation in which material also flows back onto 
the galaxy along the disk plane. The top right panel of figure 
[TJshows a cross section of the Nhi column density distribution 
m such a model viewed edge on with ^ = 45°. Wobs(t^) is ob- 
tained as before using equation [T] We refer to these models as 
(softened). 

Figure [2] shows the averaged radial profiles for four sharp 
and four softened models (averaged over both inclination and 
azimu th) compared with th e measured radial profile taken 
from ( Bordoloi et al.|2011| ) (blue data points). By construc- 
tion, all the models have an identical average radial profile 
(solid lines). Not surprisingly, the radial profile alone gives 
no information on the distribution of the absorbing mate- 
rial. 



2.2. 2-D disk model 

We also define a 2D model consisting of an infinitesimally 
thin disk, characterized as 



NperpiR) = OL 



R 

Rrui 



exp 



R 



Rr 



; along the disk (5) 



where Nperp is column density perpendicular to the disk. 
The cut-off radius {Rcut = 37.5 kpc) is identical to that of the 
3-D models, but because of the reduced dimensionality, the 
exponent of the power-law is reduced by one. The normal- 
ization term a is again chosen similar to that done in the case 
of the 3-D models. For the 2-D disk case, the observed col- 
umn density N{b) can be measured directly without integrat- 
ing along the line of sight but with a correction for non face-on 
disks because of the increase path through the disk. The radial 
position where a line of sight intersects the plane of the disk is 
given in terms of /?, / and (/) and so, accounting for the longer 
path length through a inclined disk, we get 



sin 



cose 



cos? 



l/2^ 



X sec? 



(6) 



N{b) is converted to Wobsib) using equation The model 
is again normalized to have (Wobs) — 0.7 withinZO kpc < b < 
50 kpc. The bottom panel of figure [T] shows the projected 
Nhi column density distribution of this model viewed at an 
inclined angle. In figure [2] we over plot the average radial 
profile for the disk modelT which is indistinguishable from 
that of the 3-D models. As noted above, the average radial 
profile alone contains no information to distinguish between 
2D and 3D geometries. 

3. OBSERVABLE QUANTITIES 

In this section we describe the basic observable quantities 
for the different model geometries and diagnose which ob- 
servables can most efficiently discriminate between them. 

3.1. Nniii^ (j)) column density distribution for different models 

The first observable we focus on is the distribution of ab- 
sorption column density as a function of / and at a fixed 
impact parameter. It can be used to predict the distribution of 
Mg Il-selected galaxies in (/, (j)) once a threshold is applied in 
equivalent width and after accounting for the sin / weighting 
in random orientations. This is done in section |5l 

For the purpose of presentation, we show these in Nm to 
avoid, for the time being, the question of saturation in Mg II. 
Nhi is computed using equation [2] Figure |3] shows NHiii^(l)) 
column densities for a set of sharp cones and the disk model 
at a fixed impact parameter of Z? = 45 kpc. As expected, the 
spherical model Mi go has no variation of absorber strength 
with inclination and azimuthal angles. For all the other cases 
both the sharp bipolar cone models and the disk have the 
strongest absorption for edge on systems with / approach- 
ing 90°. However for bipolar cone models, strong systems 
are preferentially aligned near the minor axis of the galax- 
ies (i.e. (/) ^ 0°) whereas for disk models this is reversed 
and the strongest absorption is aligned along the disk of the 
galaxy (i.e. (j) 90°). It can be seen that all the models, 
except the completely spherical one, have regions where no 
absorption occurs. Quasar lines of sight with strong absorp- 
tion should avoid these regions of (/, 0), while lines of sight 
without strong absorption (e.g. below 0.1 A) should cluster 
there, unless small scale dumpiness in the Mg II gas distri- 
bution permitted such lines of sight also in regions of much 
higher average equivalent width. 

3.2. Azimuthal and inclination profiles 

By collapsing one of the angular variables, we can study 
how Nhi, and/or Wots depend on the other angular variable, 
averaged over some impact parameter range. These profiles 
can be compared with the observed profi les from the stacked 
spectra analysis of |Bordoloi et al. ( 2011| ). 

To demonstrate how different model geometry is sensitive 
to / and (j), we consider the models M30, M45, M90 and M120 
as well as the 2-d disk and the sphere (Mi go) model. For 
each model Wobs and Nhi column densities are computed as 
described in the previous section. We create two observ- 
ables to compare with the stacked spectra (1) (Wobs) within 
20 < < 40 kpc, averaged over all (j), which gives the vari- 
ation of / within the b bin. (2) (Wots) within 20 < /? < 40 
kpc averaged over all /, which gives the average variation of (j) 
within the same b bin. Figure H shows the variation of (Wobs) 
as a functi on / and for these six models. The data points are 
taken from |Bordoloi et al.| ( |20TT] ). 

The bottom panel shows the azimuthal variation of (Wobs) 
where a strong variation is seen. The solid lines show the 
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Fig. 4. — The average Nhi column density distribution for different models as a function of inclination / (top panel) and azimuth (bottom panel), both 
averaged over 20 < ^ < 40 kpc. In the top panel, the i profile is averaged over all In the bottom panel, the (/> profile is averaged over all i (dashed lines) and 
just over those more highly inclined galaxies with z > 45° (solid lines). In all cases we assume that there are no saturation effects. The data points are adopted 
from the inclination and azimuthal profiles of disk galaxies in Bordoloi et al. ( 201 1) 



profile for inclined galaxies (/ > 45°). The sharp bipolar 
cone models are always stronger along the pole of the galaxy 
~ 0°) whereas the disk model exhibits strongest absorption 
near the pl ane of the disk ~ 90°). Comparing with the data 
from Bord oloi et al.| ( |2QlT] ), the disk model can be easily ruled 
out since the azimuthal dependence is of the opposite sign. It 
should be noted that the curves in figure ffl neglect the effect 
of saturation which will act to reduce the absorption where it 
is expected to be strongest (see below). It can be clearly seen 
from the top panel that the inclination dependence is rather 
a weak diagnostic of different model geometries and most of 
the models have essentially the same flat variation with incli- 
nation. 

It may strike the reader as surprising that there is such a 
weak dependence on inclination for the 2-d disk model. How- 
ever, there are two competing effects. The increased sec / path 
through the disk produces stronger absorption. On the other 
hand, as shown in equation [6] above, lines of sight (except 
those that pass through the disk along the line of nodes de- 
fined by the plane of the disk and the plane of the sky, i.e. 
those with (j) = 90° in equation [6]) will intersect the disk at 
progressively larger distances as the inclination of the disk is 
increased, reducing the absorption. Which of these wins, i.e. 
the sign of the gradient of the overall {Wobs)(i) profile, de- 
pends on the radial profile of the disk, flatter profiles increase 
the absorption at high inclinations. 

This dependence on radial profile is explored further in Fig- 
ure [5] All the models in figure |4] have the same value of Rcut- 
We now investigate how the inclination and (j) dependence 



vary when Rcut is changed. Figure [5] shows the variation in 
column density as a function of both / and (j) for the same base- 
line models if Rcut is varied from 25 kpc to 45 kpc. Clearly 
the dependence remains much the same and each model can 
be clearly discriminated. The inclination dependence is much 
more sensitive to the selection of Rcut rendering the discrim- 
ination between different models impossible. This plot high- 
lights the fact that azimuthal dependence is a much stronger 
discriminant for Mg II absorber geometry. 

3.2.1. Effects of line satu ration 

All the simulations performed above assume that there is no 
line saturation and that we are dealing with the linear regime 
of the curve of growth. This is unlikely to be valid, as dis- 
cussed in Section 2. Figure |6] shows the inclination and az- 
imuthal dependence for different sets of models after account- 
ing for line saturation using the scheme described in Section 
2. 

The left panel shows the fiducial models M30, M45, M90, 
M120, sphere and disk as described above. We also introduce 
at this point, two hybrid models Hybl (M90 + sphere) and 
Hyb2 (M90+ Disk). 

The Hybl (M^ + sphere) model is constructed by linear 
combination of a sharp bipolar cone model with opening an- 
gle x° and the sphere model. The relative strength of the two 
components is conveniently parametrized as given below. 



n = 



Psphe, 



(7) 
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Fig. 5 . — The average Nm column density distribution for M45 (solid lines), 
M90 (dashed lines), Mi 20 (dotted lines) and pure Disk (dot-dashed lines) mod- 
els which have different radial profiles obtained by setting Rcut in the range 
25 < Rcut < 45 kpc, computed as a function of inclination i (top panel) and 
azimuthal angle ^ (lower panel). It is noticeable how the sign of the incli- 
nation dependence depends on the steepness of the radial profile, because of 
the two competing effects: the increased path length sec/ vs. the effect of 
sampling greater distances in the disk - see text for discussion. The radial 
profile has little effect on the azimuthal dependence. 

where pcone gives the density of neutral hydrogen atoms in 
the cone region and psphere gives the density of Mg II atoms in 
the rest of the model. Effectively n gives the over density of 
neutral hydrogen atoms in the cone region with respect to the 
rest of the model. ^ is a free parameter, which can be varied to 
make the hybrid models have a wide range of density contrast 
between the cone region and the rest of the model. If ^ >> 1, 
the cone region dominates the final contribution to the hybrid 
model. When ^ = 1, the density of atoms in the cone region 
and the rest of the model is the same and hence irrespective 
of the opening angle, the hybrid model has a uniform density 
distribution of a spherical model. If ^ < 1, the number density 
of atoms in the cone region is less than that in the rest of the 
model. When n < 1, the over density of neutral hydrogen 
atoms will be aligned along the major axis of the disk. As 
the opening angle increases, the hybrid model approaches a 
spherical geometry. The hybrid model is normalized to have 
the same radial profile as shown in figure |3] 

The second hybrid model (Hyb2) is constructed also as a 
combination of a sharp bipolar cone with opening angle x 
and the disk model. The two components of the hybrid model 
are weighted as defined below 



f ^cone ~^ ^disk 

1^7 



(8) 



Where Whyb is the contribution to the final hybrid model, 
Wcone is the sharp bipolar cone model with a given opening 
angle and Wdisk is the disk model. / is a free parameter which 
controls the fractional contribution of each of the cone or the 
sphere/disk component. This hybrid model also has the same 
radial profile as in figure |3] 



Physically, the two hybrid models could represent a sce- 
nario where the absorbers originate due to a combination of 
a bipolar outflow component combined with either a spheri- 
cal halo component (Hybl) or an extended disk component 
(Hyb2). The inclination and radial profiles of these hybrid 
models are plotted on the left hand panel of figure^ The right 
panel of figure |6] shows the same profiles for softened bipo- 
lar cone models M30 (softened), M45 (softened), M9o(softened) 
and Mi2o(softened). 

The upper panels of figure |6] show the inclination profiles. 
As noted above, these show very little variation of absorption 
strength with inclination /. The bottom panels show the az- 
imuthal profiles. On the left hand set of models, we show 
only those objects with / > 45° since this is where the az- 
imuthal effects are strongest. On the right hand panel (with 
fewer models) we show both the azimuthal profiles for these 
inclined objects with / > 45° (solid lines) and for all inclina- 
tions (dashed lines). The hybrid models (left panel) and the 
softened model (right panel) both well represent the observa- 
tions. The optimal opening angle and the optimal n and / for 
these set of models will be determined later in Section 6. 

4. APPLICATION TO QUASAR ABSORPTION LINE SYSTEMS 

In this section we discuss how the model predictions de- 
scribed above can be translated into observables for a study 
done with quasar absorption line systems. We show how dif- 
ferent models predict the covering fraction distribution for 
strong (Wr(2796) > 0.3 A) absorbers, and study the distribu- 
tion of azimuthal angles {(j)) of such absorbers in the data from 



Kacprzak et al.| ( 2Qllb| ) and compare them to the expected 
distribution of different models. 

4.1. Covering Fraction Distribution 

Extended Mg II absorbing gas is seen out to large impact 
parameters (Z? > 100 kpc) both from quasar absorption line 
studies and stacked spectra analysis. However both the ab- 
sorption strength and gas covering fraction declines toward 
larger radii. Here, we first examine how the gas covering frac- 
tion f(Wobs > 0.6 A) varies with b and Wobs- 

In the context of our models, the gas covering fraction at an 
impact parameter b is defined as 



/(W,,,>0.6)(/7): 



6)(b) 



Ntotib) 



(9) 



^(Wobs>o.6)(b) is the number of lines of sight that have an 
absorber with equivalent width Wots > 0.6 A at an impact pa- 
rameter b and Ntot is the total number of lines of sight at the 
same impact parameter. 

Figure |7] shows the covering fraction distribution of Wobs ^ 
0.6 A absorbers as a function of b for models M30, M45, M90, 
M120, Sphere, Disk and M45 (softened). This plot is similar 
to a covering fraction distribution of Wr(2196) > 0.3 A for 
a quasar absorption line study. But as the doublet ratio of 
Wr(2796) ~ 0.3 A systems might be between 1 and 2, there 
is indeed some leeway in this conversion and comparison. As 
expected the covering fraction for the Sphere model is a step 
function with unit covering fraction below and zero cover- 
ing fraction above b = 32 kpc. Although all the sharp bipolar 
cones have exactly the same average radial profile of total ab- 
sorption, the covering fractions exhibit a range of radial pro- 
files and approach the spherical case with increasing opening 
angle. While the disk and the softened bipolar cone models 



Modeling the distribution of Mg II absorbers around galaxies 



9 



- % — \ — % -^120 % Disk-Hybl:(M^o + %o) -Hyb2:(M^o + Disk) 



M^^(softened) M^j(softened) M^p(softened) M^^q^^^''^^"^^) 



2.5 

2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
- Face-on 20kpc<b<40kpc Edge-on - 



2.5n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T 

Edge-on - 
Inclination - 




Face-on 20kpc<b<40kpc 



40 50 
Angles 



40 50 
Angles 



Fig. 6. — The mean (Wots) distribution for different models within 20 < < 40 kpc including the effects of the modelled line saturation. The left panels show 
the sharp bipolar cone models, the sphere, the disk and the hybrid models that are combinations of the Mgo and the disk or the sphere. The right panels show the 
softened bipolar cone models. The upper panels show the inclination profile, averaged over all azimuths. The lower panels are inclination-averaged azimuthal 
profiles. The (solid line) azimuthal profiles in both the left and right lower panels are averaged over i> 45°. The dashed lines in the bottom right panel are 
averaged over all inclinations. The data points are adopted from the inclination and azimuthal profiles of disk galaxies in jBordoloi et al.,^2011y 



exhibit different behavior in detail, it is hard to distinguish 
between these different geometries purely based on their cov- 
ering fraction as a function of impact parameter. For all the 
models, f(W > 0.6) ^ as /? > 80 kpc. Observations indi- 
cate that the covering fraction does not go to zero as /? > 80 
kpc, for quasar absorption line studies ( [Nielsen et al. 2012). 
But the models studied here are for isolated systems and at 
higher impact parameters contribution from filamentary ac- 
cretion, satellite galaxies etc will become significant. 

Clearly, the radial dependence of covering fraction is not 
very sensitive to model geometry. The variation of gas cov- 
ering fractions with inclination and azimuthal angles offer 
greater discriminatory power. Figure[8]shows the gas covering 
fraction for Wots > 0.6 A absorbers as a function of inclination 
angle (/) and azimuthal angle (0) within 20 < /? < 50 kpc. The 
(/) dependence of covering fraction is quite obvious from this 
figure. The sharp bipolar cone models show very high cover- 
ing fractions near the galaxy pole ((/) 0) and which gradu- 
ally decrease as the cone opening angle increases. They ap- 
proach the Sphere model as x 180°. The disk model shows 
the opposite trend and has a very high covering fraction along 
the disk axis. The inclination dependence is not as strong as 
the azimuthal dependence. Only the 2-d disk model exhibits 
higher covering fraction for the face-on galaxies (/ = 0°) as 
compared to the edge-on galaxies (/ = 90°). 

It should be stressed that any comparison of these covering 
fraction estimates with the observations should be done with 
care. These covering fraction estimates show how different 
Mg II absorption distribution geometries manifest themselves 
in terms of the observed covering fractions relative to each 
other. The models studied here represent isolated galaxies and 
the CGM is a more complex environment. These models of 
isolated systems do not include the effects of filamentary ac- 



cretion, satellite galaxies or LMC like streams or high veloc- 
ity clouds (HVCs). Further, we assume a smooth distribution 
Mg II absorption, and in reality the CGM is more patchy. At 
high impact parameters these structures will start to contribute 
significantly to the overall covering fraction. Hence the com- 
parison with the observed covering fraction should be done 
with these caveats in mind. 

4.2. Distribution of Azimuthal Angles in quasar data 
In this section we discuss the azimuthal dependence of a 
set of stron g (Wr(2196) > 0.3 A) Mg II absorbers observed 
in Kacprza k et al.| ( 2Qllb|) . [Kacprzak et al. (2011b) had a 
sample of 40 Mg II absorbing galaxies, out of which 29 have 
a rest frame equivalent width 1^^(2796) > 0.3 A. One out of 
these 29 galaxies is completely face on (/ < 5°) and it be- 
comes rather difficult to assign an accurate azimuthal angle. 
Hence we exclude this galaxy from our analysis. We shall 
focus our analysis on this sample of 28 galaxies classified as 
"strong" absorbers. Galaxy morphology, inclination and az- 
imuthal angles were obtained by using two-dimensional fits 
of a (co-centered) bulge and disk model. We refer the reader 
to Kacprzak et al.| ( |201 Ibj ) for details of methodology and data 
reduction details. 

Figure |9] shows the distribution of "strong" Mg II absorbers 
as a function of their azimuthal angle ((/)) and impact param- 
eter (b). 72% of all the "strong" absorbers (20 out of 28) 
that are within b < 120 kpc lie within 50° of the semi minor 
axis of the galaxies i.e. to the left of the dashed vertical line 
of figure [9j as compared with 55% that would be expected 
for a uniform azimuthal distribution. More significantly, for 
Z? < 38 kpc (below the horizontal dashed line) all of the nine 
"strong" absorbers lie within 50° of the galaxy semi minor 
axis. This evidence for a strong azimuthal dependence sup- 



10 



Bordoloi et al. 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




50 60 70 

Impact Parameter (kpc) 

Fig. 7. — The covering fraction for Wobs > 0-6 A absorbers as a function of impact parameter in different models, all based on idealized smooth 3-d density 
dtributions. The spherical model shows a step function behavior as expected. The behavior of the sharp bipolar cones vary with opening angle and approach 
the spherical case with increasing opening angle. The disk and the softened bipolar cone models exhibit different behavior but it is hard to distinguish between 
different models using this diagnostic. 
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Fig. 8. — Covering fraction for Wots ^ 0-6 A absorbers as a function of i and (j) for different models. Clearly the azimuthal dependence of covering fraction is 
a much stronger discriminant of absorber geometry although inclination dependence has a weak correlation with different models. 
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Fig. 9. — The distribution of "strong" Mg II absorbers ( Wr(2796) > 0.3 A) in azimuthal angle and impact parameter b (adopted from Kacprzak et al. 2011b 
72% of the absorbers lie within 50° of the semi minor axis of the galaxies, i.e. lie to the left of the vertical dashed line. The horizontal dashed line marks b < 38 
kpc. Within this distance, all of the absorbers lie within 50° of the semi-minor axis of the galaxy. 




Fig. 10. — Cumulative distribution of azimuthal angles ((/>) moving away from the projected minor axis for different 3-d models and for the quasar absorption 
lines with (Wr(2796) > 0.3 A, adopted from ( Kacprzak et al. 201 lb I. The data clearly rule out a spherical or a disk-like geometry alone. However a hybrid 
cone-i-disk or cone-i-sphere model or a softened bipolar cone structure is consistent with the data, depending on the opening angle of the cone. 
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ports the hypothesis that out to b ~ 40 kpc the "strong" Mg 
II absorbers are originating in bipolar galactic winds, as indi- 
cated by the mea surement of total Mg II absorption in |Bor- 
Idoloi et al.| (|20TT]). 

To compare this dataset with models we plot the cumulative 
distribution of (j) of the Mg Il-selected galaxies and a set of 
models in figure |T0| The gray squares represent the cumula- 
tive distribution function for the absorbers detected at /? < 38 
kpc. Clearly this dataset favors a more cone like distribution 
of absorbers. The disk and sphere models do not match the 
data. Even for the absorbers at /? > 38 kpc (black circles), 
the quasar absorption line data immediately rule out the disk 
model. Also M30 and M45 models do not match the data. The 
data can however be reproduced by the softened bipolar cone 
model or a hybrid model. At low impact parameters, the data 
are consistent with models of small opening angle whereas 
the high impact parameter data are consistent with high open- 
ing angle models. In the next section, we quantify the models 
supported by the quasar absorption data as well as the inte- 
grated spectra data. 

5. JOINT ANALYSIS OF STACKED SPECTRA AND QUASAR 
ABSORPTION LINE SYSTEMS 

In this section, we carry out a joint analysis of the two inde- 
pendent and orthogonal probes of Mg II absorbers discussed 
above, namely the stacked spectra data of Bordoloi et al. 
([2011) and th e quasar absorption line probes from Kacprzak, 
|etal.||2011b| ). 

We chose a grid of models with opening angles within 
1° < X ^ 170° for both the sharp bipolar cone and the soft- 
ened bipolar cone models. We create, as above, two sets of hy- 
brid models combining a sharp bipolar cone model first with 
a sphere model and second, with a disk model as described 
in the previous Section 3.2.1. We examine the cone+sphere 
hybrid model in the over-density (equation [7]) range of 0.1 < 
^ < 10 and the cone+disk model in the range of 0.1 </< 10 
(equation [8]). For each model the overall normalization factor 
a is assigned so that {Wots) ^ 0-7 A within 20 kpc < b < 50 
kpc. 

To compare the models with the integrated spectra data, we 
compute the mean azimuthal profile for each model within a 
given bin of impact parameter. The mean azimuthal profiles of 
the models and the measurements are compared by computing 
the value for each model, taking the minimum value as 
the best fit parameters for each particular set of models. 

The quasar absorption line data are compared with the same 
set of models in the same b bins by performing a two di- 
mensional Kolmogorov-Smirnov test ( [Pea cock 1983; FasanoJ 
|& Franceschini 1987; Press 2002|). To perform the two di- 
mensional Kolmogorov-Smirnov test we chose all "strong" 
absorbers (above 0.3 A) within a fixed b range and compare 
the distribution of inclination and (j) angles with that expected 



for different models. Figure 1 1 shows the distribution of / and 
^ for the"strong" absorbers m six models within 20 <b <50 
kpc and the d istribution of "strong" absorbers from | Kacprzak| 



et al. ( 2011b ). In this b range, the pure (sharp) bipolar cone, 
sphere and the pure disk models are all ruled out by the 2-D 
KS test at 5% significance level. The softened bipolar cone, 
the hybrid cone+sphere and the hybrid cone+disk models are 
all acceptable (cannot be ruled out) by the 2-D KS test. This 
reflects the fact that while most of the absorption is above the 
poles, some significant amount is required at large azimuth 
angles. 

To estimate the model parameters which best represent the 



observations, we compare the observations with the model 
predictions within b < 40 kpc. We first focus on the soft- 
ened bipolar cone model and compare the models within 
20 kpc < b < 40 kpc with the observations. Figure [12] shows 
the joint constraints from both the integrated spectra3ata and 
the quasar absorption line data. The red curve shows the like- 
lihood curve from the integrated spectra data and the gray 
hashed region gives the opening angles allowed by the 2-D 
KS test at 10% significance level (i.e. parameters outside of 
this region are excluded with more than 90% confidence). The 
black hashed region is the joint constraints from both these 
probes which are obtained by simply taking the product of the 
probabilities. This analysis gives an opening angle of 45^1^ 
degrees that is allowed by both the observations for a softened 
bipolar outflow model at 90% confidence. 

We now shift our focus on the two hybrid models. We first 
examine the hybrid model having a cone and a sphere compo- 
nent. As described above we compute the x^ values, by com- 
paring each model with the integrated spectra. The allowed 
confidence contours for the two parameters of interest, n and 
opening angle are shown in figure [13] The contours show the 
68%, 90% and 95% confidence intervals of the allowed model 
parameters. The blue contours are allowed regions in parame- 
ter space obtained from the quasar data at 95% and 90% con- 
fidence level. The hashed regions give the joint constraints 
on opening angle and n at 95% and 90% confidence levels 
respectively. The left panel shows the allowed model parame- 
ters within 20 kpc < b < 40 kpc. At these impact parameters, 
both the integrated spectra data and the quasar absorption line 
data favor models where there is an enhancement of Mg II 
absorption within ~ 100° of the bipolar region. The opening 
angle allowed by both the datasets within 40 kpc is 105;^^q° 
and the over-density n is 2;^q 2. 

We now look at the second hybrid model having a cone and 
a disk component. Figure [14 shows the constraints on the 
allowed model parameters, tEe plotting scheme is the same 
as described above. The close in data (left panel, 20 kpc < 
b < 40 kpc) show that both the quasar absorption line systems 
and the integrated spectra are in agreement with a cone like 
geometry with an opening angle of ^ 110°. The quasar ab- 
sorption line data allow for a range of opening angles which 
are smaller than those allowed by the integrated spectra data. 
We put joint constraint on the opening angle at 1 10° ± 5° and 
the fractional contribution to the cone component is found to 
be/=6!t. 

The Mg II absorption around galaxies show strong az- 
imuthal asymmetry for b < 40 kpc. But at higher impact pa- 
rameters, this azimuthal asymmetry is much weaker as shown 
in Bordoloi et al. (2011). We now investigate the models 
that can best represent the observations within 41 kpc < b < 
80 kpc as shown in the right panel of figure JT3] We see 
that the models with n= 1 i.e. spherical distribution of Mg 
II gas around the galaxies, cannot be ruled out by either of 
the datasets. Both the datasets are insensitive to the very low 
and high opening angles since such very sharp features pro- 
duce little effect on the overall model. By combining the two 
probes together we see that the spherically symmetric distri- 
bution with ^ ~ 1 is strongly preferred. 

Furthermore, the joint models obtained for the two impact 
parameter ranges at 90% confidence do not overlap. This is 
evidence that the distribution of Mg II absorbers changes from 
a highly asymmetric distribution at low impact parameters to 
a less asymmetric distribution, consistent with a spherical dis- 
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Fig. 11. — The distribution of "strong" Mg II absorbers ( Wr(2796) > 0.3 A) within b <50 kpc as a function of incHnation and azimuth. The panel at the 
bottom shows the data of Kacprzak et al. ( 201 lb). The distribution of strong absorbers in the models with a sharp bipolar cone (M45, top left), with a sphere (top 
right), a disk (third panel from top), a softened bipolar cone (fourth panel from top), a hybrid cone+sphere (fifth panel from top) and a hybrid cone+disk (sixth 
panel from top) are also shown. A two dimensional KS test is performed to compare the distribution of absorbers in the models with the observations. 
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tribution, at high impact parameters. 

A similar conclusion is reached from considering the sec- 
ond hybrid model (cone+disk) at the larger impact parameters 
(right panel, figure 14). The quasar data are quite degenerate 
in / and opening angle (requiring / to reduce as the opening 
angle is narrowed) but are consistent with the stacked galaxy 
analysis. For allowed small opening angles there is signifi- 
cant contribution from the disk component. As the allowed 
opening angle increases the cone component becomes more 
and more dominant. Combining these two probes, we get 
joint constraints on the opening angle = 130:^45 ° / " ^-i.i • 
This best fit model also represents an almost spherical geome- 
try with very large opening angle for the cone component. As 
before, the region allowed at small impact parameters does 
not overlap with the region allowed at larger radii. 

We conclude that the distribution of Mg II gas at low im- 
pact parameters is not the same as that found at high impact 
parameters. This radial trend is hard to see on either the in- 
tegrated spectra or the quasar absorption line data alone be- 
cause individually, either probe is not sensitive enough to rule 
out the same regions of parameter space at different b (fig- 
ure [13] [14]). However the joint analysis reveals that the Mg II 
gas IS distributed preferentially along the minor axis of disk 
galaxies with an opening angle of 100° and an overdensity 
of ^ ~ 2 - 3 and f 6-1 within /? < 40 kpc. At high impact 
parameters the distribution is no longer aligned only along the 
minor axis and a spherical distribution is also consistent with 
both the dataset. 

These models quantify the radial changes in the distribution 
of Mg II absorbers around galaxies. With the best fit models 
we can present a more clearer picture of the radial distribu- 
tion of Mg II absorbing gas. At close impact parameters (e.g. 
/? < 40 kpc), the Mg II absorption is primarily aligned along 
the minor axis of the galaxy as they are arising due to mate- 
rial entrained in star formation driven outflows. As we move 
outwards, the distribution of the absorbers are consistent with 
either a spherical distribution or a cone with a large open- 
ing angle augmented by an extended disk component which 
might be the signature of in falling material, perhaps as part of 
a galactic fountain or a contribution from material falling in 
from larger distances. The more spherical distribution could 
also reflect the contribution of satellite galaxies. 

6. CONCLUSION AND DISCUSSIONS 

In this paper we address the principle question of the distri- 
bution of cool Mg II gas observed as absorption line systems 
in the spectra of background sources. We present a suite of 
3 -dimensional models of Mg II absorption to study the az- 
imuthal and inclination dependence of observational quanti- 
ties, including the integrated column density and the proba- 
bility of intersecting a region above a given column density 
threshold. We compare these model distributions to observa- 
tions of Mg II absorption that we have previou sly obtained by 
stacking the spectra of background galaxies (fBordoloi et al. 
201 1| ) and by analyzing, in this paper, the distribution Mg II 
absorbers in quasar spectra as a function of the impact param- 
eter and relative galaxy orientations ( Kacprzak et al. 2011b). 
Both these orthogonal probes of cool ionized Mg II gas, allow 
us to put joint constraints on the distribution of Mgll gas. The 
main findings of this analysis can be summarized as follows. 

1 . The azimuthal dependence of the observable quantities 
is a much stronger diagnostic of the 3-d geometry than 
the dependence on inclination. 



2. A composite model that consists of a simple bipo- 
lar cone plus a spherical or disk component or a sin- 
gle softened bipolar distribution can well represent the 



azimuthal dependencies observed in ,Bordoloi et al. 
( |2Qn| at < 40 kpc. 



3. The stro ng Mg II absorbers observed in [Kacprzak et al.| 
([201 IW are also asymmetrically distributed. Overall, 
72% of these the absorbers are within 50° of the galaxy 
semi minor axis. Within Z? < 38 kpc, "all" the strong 
absorbers are within 50° of the galaxy semi minor axis. 
Comparing these absorbers with models and plotting 
the cumulative distribution of the azimuthal angles, we 
see that the data is consistent with the distribution of az- 
imuthal angles of bi-polar models and any purely spher- 
ical or disk model can ruled out at low impact parame- 
ters. 

4. A joint analysis combining the integrated spectra data 
and the quasar absorption line data shows that within 
40 kpc both these datasets can be well reproduced by 
a softened bipolar cone model of opening angle 45;^9^ 
degree or hybrid models consisting of a cone and a 
sphere component with n is and an opening an- 
gle of 105;^^Q degree. A second hybrid model with a 
cone and a disk component is also consistent with the 
observed data with an opening angle of 110° ± 5° and 

5. These analyses show that the Mg II absorbers at low 
impact parameters are primarily distributed along the 
minor axis of disk galaxies and are originating in galac- 
tic winds, with only a small fraction might be aligned 
along the major disk axis. In rough terms, the number 
densities of Mgll above the poles of the galaxies are 2-3 
times higher than in the planes of the disks. 

6. At larger impact parameters (41 kpc < /? < 80 kpc), the 
distribution of Mg II gas is more symmetric than that 
observed at smaller impact parameters. A joint analysis 
of the integrated spectra and the quasar absorption line 
data shows that the model best representing both the 
datasets is a more or less spherical distribution of Mg II 
gas. 

7. The clear change in the distribution of Mg II absorption 
with increasing impact parameter suggests that the ori- 
gin of the absorption may also change as we move out 
in radius from the galaxies. The more symmetric distri- 
bution at large radii may reflect an increased contribu- 
tion from in-falling material, or even satellite galaxies. 

We also present predictions for the expected covering frac- 
tion of strong absorbers as a function of impact parameters 
for different absorber geometries. The covering fraction as a 
function of impact parameter is not a strong discriminant of 
absorber geometry. 

We predict the expected inclination and azimuthal depen- 
dencies for individual QSO sight lines within 20 kpc < /? < 50 
kpc for a set of given models. The covering fraction of Mg 
II absorbers are most sensitive to galaxy azimuthal angles. 
Amongst the covering fraction profiles, the azimuthal profile 
of gas covering fraction is the most sensitive probe to discrim- 
inate different absorber geometries. 
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Opening Angle (Degree) 

Fig. 12. — The joint constraints on the opening angle of the softened bipolar cone model at 20 < < 40 kpc from both the stacked spectra analysis and the 
quasar absorption line studies. The red solid line is the normalized likelihood curve obtained from the minimization of the model predictions to the azimuthal 
profile in the stacked spectra analysis. The maximum likelihood opening angle of about 45° is consistent with the opening angles allowed by performing a 2D KS 
test of the model angle distributions with the quasar absorption line data (shown as the open hashed region). The close hatched region shows the joint constraint 
at 90% confidence level. 
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Fig. 13. — Joint constraints on the hybrid Cone+Sphere model from both stacked spectra analysis and the quasar absorption line studies on two parameters 
describing this model - these are n, the density contrast above the pole and in the plane of the disk, and the cone opening angle. The 68%, 90% and 95% 
confidence contours obtained from fitting of the model predictions to the azimuthal profile in the stacked spectra analysis (gray contours). The blue contours 
show the parameter space allowed by the quasar absorption line data at 90% and 95% confidence level, obtained from a 2D KS test analysis. The hashed region 
on both panels show the parameter space allowed by both datasets at 90% and 95% confidence level. The left panel is for absorption within 20 ^pc < ^ < 40 kpc, 
and the right panel is for absorption within 40 kpc < b < ^0 kpc. The differences show the pronounced weakening of the azimuthal asymmetries with radius. 
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